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Abstract 

In this paper we study double scalar neutrino inflation in the minimal supersymmetric seesaw 
model in light of WMAP. Inflation in this model is firstly driven by the heavier sneutrino field N2 
and then the lighter field A'"i. we will show that with the mass ratio 6 ^ M2/M1 ^ 10 the model 
0> \ predicts a suppressed primordial scalar spectrum around the largest scales and the predicted CMB 

TT quadrupole is much better suppressed than the single sneutrino model. So this model is more 
. favored than the single sneutrino inflation model. We then consider the implications of the model 

\ on the reheating temperature, leptogenesis and lepton flavor violation. Our results show that 

I the seesaw parameters are constrained strongly by the reheating temperature, together with the 

m ■ 

■ requirement by a successful inflation. The mixing between the first generation and the other two 



generations in the right-handed neutrino sector is tiny. The rates of lepton flavor violating processes 



Q^' in our scenario depend on only 4 unknown seesaw parameters through a 'reduced' seesaw formula, 

' besides C/e3 and the supersymmetric parameters. We find that the branching ratio of — > 67 is 



generally near the present experimental limit, while Br(r — > ^7) is around 0(10 — 10 '^^ 
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I. INTRODUCTION 



It is widely accepted today that the early universe has experienced an era of accelerated 
expansion known as inflation Inflationary universe has solved many problems of the 
standard hot big-bang cosmology, such as the flatness and horizon problems. In addition, 
it provides a causal interpretation for the origin of the density fluctuations in the Cosmic 
Microwave Background (CMB) and large scale structure (LSS). 

nn 

Among current inflation models, sneutrino chaotic inflation [2, |3[ is one of the promising 
physical candidates where inflation is driven by the superpartner of the right-handed (RH) 
neutrino. In this scenario, no extra inflaton scalar field is needed, besides the RH sneutrinos, 
which are _ to exp.a,„ *e ti.y ,„a.Q i„ the .ini.a. .upe.sy,n.et.c 

seesaw mechanism j5|. Baryon number asymmetry via leptog also be easily 

realized in this framework. 

The single sneutrino inflation model predicts a near scale invariant primordial power 
spectrum. Despite the fact that the scale invariant primordial spectrum is consistent with 
current Wilkinson Microwave Anisotropy Probe (WMAP) observations 7], it is noted that 



there might be possible discrepancies between predictions anc 
and smallest scales. WMAP data show a low TT quadrupole 



observations on the largest 
as previously detected by 



COBeIqI. In Ref.Q Peiris et al. find that WMAP data alone favor a large running of the 
spectral index from blue to red at ^ 1.5a with dns/dlnk = — 0.077io;Q52. When adding LSS 



data of 2DFGRS^^ the running is more favored with dns/d\nk = —0.075 



+0.044 
-0.045- 



The most proper way to get the shape of the spectrum from observations should be 



12 



141. A detailed reconstruction of the power 

n 

12| shows that a running of the index is favored. Ref. [13|| 



the primordial spectrum reconstruction 
spectrum by Mukherjee and Wang 
reconstructs the primordial spectrum with WMAP data and the shape of the matter power 
spectrum from 2DFGRS^|. The authors attribute the need for the running to the first 
three CMB multipoles / = 2,3,4. They introduce power-law spectrum with a cut at large 
scales and find a non- vanishing cutoff is favored at ^ 1.5cr. 

The statistical level of the low CMB multipoles has been discussed wide^ 




many models have been built to achieve the suppressed CMB multipoles 
though the confidence level of spectral index run ning is not very high, if stands, it would 
severely constrain inflation model buildings jiol. l2ol. I21I and the single field sneutrino chaotic 
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inflation model would be in great challenge^. 

Recently we have considered a double inflation model 



m 



V{<P,,ct>2) = \ml4>l + \ml<Pi , (1) 

where inflation is driven flrstly by the heavier inflaton <^2; then the lighter fleld 0i. But 
there is no interruption in between. This model solves the problems of flatness etc. and 
generates a primordial spectrum suppressed at certain small k values. The CMB quadrupole 
predicted can be much lower than the standard power-law ACDM model. Recently, it is 
shown by Kamionkowski et al.'2^] that the cross-correlation between the CMB and an all- 
sky cosmic-shear map will be enhanced by such a primordial spectrum, and this may be 



observable at 2 — Sex 
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. The suppressed CMB multipoles can also lead to many other 
3es3- 



observable consequences 

In the present work, we consider the case that the two inflaton flelds consist of the two 
lighter sneutrinos, A'^i and N2 in the minimal supersymmetric seesaw model, while the heav- 
iest one, A'^s, does not contribute to inflation. By fitting the resulted primordial spectrum 
to the WMAP data in the next section, we get the preferred two sneutrino masses. Mi and 
M2. We find that the double sneutrino model is more favored than the single sneutrino 
model at about 1.5cr level. In section III, we first present our parameterization of the seesaw 
model and then analyze the implications of this model on the reheating temperature, lep- 
togenesis and lepton fiavor violation, etc. We find the reheating temperature, constrained 



by the gravitino problem 27| to be below 0{lQ^'^GeV), gives very strong constraint on the 
seesaw parameter space and our analysis is greatly simplified then. Different from a random 
sampling on the 9-dimensional unknown seesaw parameter space in Ref. we can show 
the seesaw parameter dependence of the predicted lepton flavor violating rate explicitly. 
Our analysis shows that there is no direct connection between leptogenesis and LFV in this 
model. Non-thermal leptogenesis is easily to be achieved via the sneutrino inflaton decay. 
Only hierarchical neutrino mass spectrum at low energy can be produced and the neutrino- 



less double beta decay 28| can not be explained by the effective Majorana neutrino mass in 
the model. 



^ Several authors in the hterature have fitted WMAP using different codes or adding various CMB and LSS 
data, they give consistent results^, Q but with less hints for running of the spectral index. 
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II. DOUBLE CHAOTIC SNEUTRINO INFLATION 



The evolution of the background fields for double sneutrino infiation is described by the 
Klein-Gordon equation^: 

0, + 3iJ0, + V^^ = Q, (2) 

and the Friedmann equation: 



H-^ - ' — 



(3) 



a 3 

where J = 1,2, a is the scale factor, the dot stands for time derivative and Vx = dV/dx. 
Defining the adiabatic field a and its perturbation as js^: 

& = (cos6')0i + (sin^)^2 , 
5a = {cos9)S(f)i + {sin 9)6 (f)2, (4) 



with 



cos6' = , , sint' = , . (5) 

/ . 9 To ' / . 9 r9 \ / 



2 -2 / ■ 2 -2 

>1 + 02 V 01 + 02 



The background equations and (jH)) become 

a- + 3Ha + K = , (6) 

where = (cos6')V^j + (sin6')V02- We assume an adiabatic initial condition between the 
perturbations 501 and 502: 

501 ^ 502 

^ 01 02 

As shown in Ref.|22i], if the initial perturbation is adiabatic, it will remain adiabatic on large 



scales during infiation. In this sense, infiation is equivalently driven by a single infiaton a 
with the effective potential V{(t) = V^(0i) + V^(02)- The basic picture of inflation and 
perturbation in our model is: the heavy inflaton 02 rolls slowly down its potential and starts 
to oscillate when the Hubble expansion rate is around its mass H ~ M2, while 0i remains 
slow rolling and ^(0i) comes to dominate the inflaton energy density. Hence, inflation is 
not suspended during the transition. 



^ To be consistent with the usual convention, in this section, we use 0/ to represent the inflatons, the 
sneutrinos here, instead of the symbol Nj. 
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FIG. 1: Effective potentials V{a) together with V{(f)i). The horizontal axis is the value of inflaton 
01 or a, in unit of Mpi. The vertical axis delineates the inflaton potential, in unit of M^^. 

The effective potential V{a), as well as the background evolution, is determined by the 
initial values of 0i, 02 (i-e. (pu and 02j) and their masses Mi and M2 (or equivalently Mi 
and r = M2/M1 ). As the heavier inflaton oscillates, | 02 \oc a~i, V^(02) oc a~^, and becomes 
negligible, one has a = 0i and V{a) = V{(f>i). Therefore, the value of a can be set equal to 
01 and they would have the same potentials. In FigQwe show the effective potential V{a) 
as well as V^(0i). V{a) becomes sharper as r increases and the initial value of 0i would also 
change the shape of the effective potential. 

We notice, from Fig. ^ a achieves a large value during the transition time and the scalar 
power perturbation is suppressed via the slow-rolling(SR) formula Ps oc (^)^- The SR 
parameters e and 6 during the transition are 

e.-|=4.G(|f4^. (8) 
5 
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V(a), r=12 




and 

... ... .2 

^ _ 0101 + 0202 _ 302 , , 



i7(0r + 02') 0r + 



^2 



■ 2 _Q 

We notice that when 02 oscillates, p^^ ~ 02 oc a , e and —5 reach their local maximum 
values. One can also find the maximum value {—6)max > ^max- In the extreme limit when 
V{(f)i) is negligible during the transition one has {—S)max = 3 and e^ax = 1-5. Regarding the 
fore-mentioned four parameters, the ratio r = M2/M1 and the initial value of 0i, determine 
the locations and values of —Smax and emax- The maximal values are mainly determined by 
r. If the ratio r is too small (e.g. 1 < r ^ 3) the above picture cannot be realized because 
both fields would take effect during inflation and neither is negligible. While r is too large 
(e.g. r > 11) one gets 1 + e + 6 < during the transition and superhorizon effectsj^^ would 
take place. The perturbations do get suppressed at some smaller k but enhanced around 
certain larger k values. Under such circumstances the whole effect might be negative to 
achieve small CMB TT quadrupole. The need that Ps{k) be suppressed at small k requires 
some tuning of the initial value of 0i. Mi determines the amplitude of the perturbation and 
is normalized by the current observations. The initial value of 02 is arbitrary with a weak 
prior to provide enough number of e — folding to solve the flatness problem. 

As our model parameters lie in the region where SR approximation does not work well, 
we calculate the primordial scalar and tensor spectra using mode by mode integrations jisl. 



^,|31|. We denote the scale where Ps arrives around its local maximum as kj and tune 



the initial 0i to get N{kf) ~ 55. In Fig. |21we show the numerical results of the scalar and 
tensor spectra for r = 3.5, 9 and 11.5. One can see that, for r = 3.5, the spectra is almost 
featureless while well suppressed scalar spectra have been generated for r = 9 and 11.5. For 
the example of r = 11.5 is enhanced around kj due to the superhorizon contributions js^l. 
We then flt the resulting primordial spectra to the current WMAP TT and TE data. As 



shown in Refs. 



n models one cannot know the exact values oi kj due 
to the uncertainty in the details of reheating. So In/cj is another parameter in our model. 
Our fltting is simflar to Ref.H: We flx Qbh^ = 0.022, fi^/i^ = 0.135, = 0.17, Qtot = 1 
[15I 1 and set Qa and In A;/ as free parameters in our flt. Denoting kc = 7.0 x 70./3/10^ ~ 
1.6 X 10^^ Mpc~^, we vary grid points with ranges [0.68,0.77], and [—3,5.] for and 
\n{kf/kc), respectively. M2/M1 varies from 3.5 to 12 in step of 0.5. At each point in the 
grid we use subroutines derived from those made available by the WMAP team to evaluate 
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FIG. 2: Primordial scalar Pg and tensor spectra Pg for r = 3.5, 9 and 11.5. The overall amplitude 
can be normalized by WMAP. 

the hkehhood with respect to the WMAP TT and TE data 34]. The overall amplitude of 
the primordial perturbations has been used as a continuous parameter. 

In Fig. 121 we plot the resulting values as functions of r and ln{kf/kc). The contours 
shown are for Ax^ values giving 1.1, 2, and 3 a contours for two parameter Gaussian 
distributions. As the location is rather hard to be fixed at exactly N{kf ) = 55, the figure is 
not very smooth as expected. Our main intention is to see how the primordial spectrum with 
a feature is favored by WMAP. This can be also seen in the one-dimensional marginalized 
distribution of ln{kf/kc) for each r. To see clearly how the feature is favored, we do not 
marginalize over r and show some of them in Fig. HJ For r = 3.5, kf ^ is favored and 
when r = 7, ln{kf/kc) = 3 is favored at around 2a. /c/ ~ is excluded at less than la for 
r = 4.5 where Ps is not suppressed enough around kf. While for r ^ 11, Ps is enhanced 
around kf and although nonzero kf is favored for shown examples, fcj ~ is excluded at less 
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FIG. 3: Two-dimensional contours in the r- ln{kf /kc) plane for our grids of model, kc ~ 1-6 x 10 ^ 
Mpc~^ The regions of different color show l.lo", 2 and 3a confidence respectively. 

than la. We find that, for 6 ^ r < 10, nonzero kf is favored at ^ 1.5a. For the investigated 
parameter space with 3.5 ^ r ^ 12 we have Ps{0.05/Mpc) = 2.46 ~ 2.59 x 10~^ at 2a leveL 
This gives Mi ~ 1.7 x lO^^ GeV. 

A detailed analysis give, .he e-fo.ds nu.be. NW> before e.d of ..fla.ionHQ: 



iV(^) = 60.56 -ln/.-ln-^-ln-^ + ln^--ln^ 



(10) 



Pend PrH 

where p{k), pend denote the inflaton potential at = aH and at the end of inflation re- 
spectively, pnH is the energy density when reheating ends, resuming a standard big bang 
evolution. Since in our case there is a preferred scale ln{kf/kc) while N{kf) is fixed around 
55, the reheating energy may be determined by the current observations. However, one can 
see that the location of fc/ is mainly determined by the initial value of 0i. Once the initial 
01 changes, N{kf) will change and the resulting pfin would be different. We show the case 
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FIG. 4: One-dimensional marginalized distributions of \n{kf/kc) for r = 3.5, 4.5, 6.5, 7, 9.5, 11 
and 12. 

in Fig. El as an example. For r = 8, (pu = 3.2 and 3.3Mpi lead to N{kf) = 54.34 and 59.06 
respectively. We get Ps ~ 2.5 x 10"^ at 0.05 Mpc'S the resulting \n{kf/h) = (-1.7, 1.3) 
and (0.1, 1.9) at 2a respectively. We also have h ^ 0.73, p3(0.05/Mpc) ~ 1.8 x lO^^Mp^ and 
pirf ~ 5.1 X IQ-^Mpi. Taking these to the models we get pIh = (8.5 x lO^GeV, 6.9 x lO^GeV) 
and (2.6 x lO^^GeV, 5.8 x lO^^GeV) at 2a for the two different (pu. Therefore, the reheating 
temperature is fully correlated with initial 0i in this model. 

We get our minimum = 1429.1 when r = 8.5 and ln{kf/kc) = 2.4. When compared 
with the standard power-law ACDM model, we have minimum = 1432.7 and Ax^ = —3.6. 

For the single field chaotic inflation we get minimum = 1432.9, with Ax^ = —3.8. 

1 

However, in the sneutrino inflation, we have to set ^ 10^'^ GeV due to the gravitino 
problem[27[. In this case, we get N{kc) ^ 55.5 and minimum = 1433.2, which gives 
Ay^ = —4.1. In addition, there are only two parameters, the mass and hi{kf/kc), in the 
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FIG. 5: One dimensional likelihoods of r=8, (jju = 3.2 and 3.3 Mpi. 

single field sneutrino inflation model. This indicates our double sneutrino inflation is favored 
at ~ 1.5(T by WMAP than the single field sneutrino inflation. In Fig. IHlwe show the resulting 
CMB TT multipoles and two-point temperature correlation function for single and double 
field sneutrino inflation in our parameter space. One can see that the resulting CMB TT 
quadrupole and the correlation function at ^ 60° are much better suppressed in the double 
sneutrino inflation than in the single sneutrino model. In fact, the spectrum of the single 
field sneutrino inflation is equivalent to that in our double case with r = 1 and (pu = (f)2i- 
In this sense we get 6 ^ r < 10 is favored at ^ 1.5o" (Ax^ ~ —2.3) than r = 1 in double 
sneutrino inflation. 

Finally, it is worth mentioning that we have also considered a double inflaton model with 
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FIG. 6: CMB anisotropy and two-point temperature correlation function for single and double 
field sneutrino inflation. Left: From left top to bottom, the lines stand for single sneutrino inflation, 
double sneutrino inflation with ln(kf/kc) = 3.0, 3.2, 3.4 and 3.6. r is fixed at 8.5. Right: From 
right top to bottom, the lines stand for single sneutrino inflation, double sneutrino inflation with 
ln{kf/kc) = 3.0, 3.2, 3.4 and 3.6 and the WMAP released data. 

quartic potential^ 

^01,02) =Ai0t + A20^ . (11) 

As we known, the quartic potential X(j)'^ is disfavored by the current WMAP and LSS ob- 
servations, because it has a larger tensor perturbation. Peiris et al. [10] fix the number of 
e-folding at 50 and find Xcj)^ infiation model is excluded at more than 3cr by WMAP and 
2DFGRS data. WMAP alone excludes infiation at more than 99% confidence level 

^ The quartic term of sneutrino is absent in the minimal supersymmetric seesaw mechanism. These terms 
can arise if the RH neutrino Majorana mass is produced in the superpotential X^NN , with $ another 
superfield whose vacuum expectation value generates the Majorana mass. 
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when ~ 50. The discrepancy between the theoretical predictions and observations comes 
mainly from the contributions of small CMB multipoles. In the double inflaton quartic 
model, the CMB quadruples can also be well suppressed and the model is also favored by 
WMAP. We fix N{kf) = 50 and run two codes, one with A2/A1 = 6400 and the other with 
= 3600 and fit the primordial scalar and tensor spectra to WMAP TT and TE data. 
We get minimum = 1427.9 and 1428 respectively. They work better than the double 
quadratic sneutrino infiation. Reheating temperature in this case cannot be restricted from 
WMAP, as shown in Ref. 3- 



III. PHENOMENOLOGY 



In the minimal seesaw mechanism, the right-handed sector is least known. However, in 
the double sneutrino infiation model, two neutrino masses Mi and M2 are constrained by 
the WMAP as shown in the previous section. In the following, we will study the phenomeno- 
logical implications of this model, including the reheating temperature, leptogenesis, lepton 
fiavor violation and neutrinoless double beta decay. 



A. Parameterization of the minimal seesaw model 



In this subsection we present our convention and parameterization of the minimal super- 
symmetric seesaw model. At the energy scales above the RH neutrino masses, the superpo- 
tential of the lepton sector is given by 

W = Yl'*H^LiE, + y;^*H2UNj + ]^MfN,N, + ^iHiH-, , (12) 

where Yl and Y^ are the charged lepton and neutrino Yukawa coupling matrices, respec- 
tively, Mr is the Majorana mass matrix for the right-handed neutrinos, with i and j being 
the generation indices. 

Generally, Yi and Y^ can not be diagonalized simultaneously. This mismatch leads to 
the lepton fiavor violating (LFV) interactions. The three matrices Yl, Yn and Mr can be 
diagonalized by 

Yl = uIYlUr , (13) 
Y^ = vIy^Vr , (14) 
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Ml = XV^MrVrX^ , 



(15) 



respectively, where U , Vl^r and X are all unitary matrices. 

We can define the lepton fiavor mixing matrix V, the analog to the Kobayashi-Maskawa 
matrix Vkm in the quark sector, as 



V = uIvl 



(16) 



V is determined by the left-handed mixing of the Yukawa coupling matrices Yl and Fjv, and 
only exists above the energy scales Mr. We will see below that this matrix determines the 
LFV effects in the supersymmetric seesaw model at low energies. 

We then rotate the bases ol L, E and N to make both Yl and Mr diagonal. On this 

basis, Yat can be written in a general form as 



Ym = VY'X 



(17) 



By adjusting the phases of the superfields, y is a CKM-like mixing matrix with one physical 
CP phase, and X has the form 

U \ (l ^ 



X = 



V 



X 



(18) 



where a, /3, p and uj are Majorana phases and X is a CKM-like mixing matrix with another 
Dirac CP phase. It is then easy to count that there are 18 parameters to parametrize the 
minimal seesaw mechanism, which include 6 Yukawa coupling constants (or mass) eigenval- 
ues in y^v and Af/j, 6 mixing angles and 6 CP phases in V and X. 

At low energies, the heavy RH neutrinos are integrated out and the Majorana mass 
matrix for the left-handed neutrinos is given by 



rriy = —rriN 



M, 



T 



R 



(19) 



where rriN = Fat?; sin /j is the neutrino Dirac mass matrix, with v being the vacuum expec- 
tation value (VEV) of the Higgs boson, rrii, can be diagonalized by 



m,. 



(20) 
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where = Uv ■ diag(l, e*^, e*^) is the MNS mixing matrix 36j| . with r], ^ being low energy 
Majorana CP phases. Ui, describes the neutrino mixing at low energies, which is different 
from the high energy mixing matrix V defined in Eq. fllfij) . From Eq. (fTIHl we can see that 
nil, is related to all the 18 seesaw parameters. However, measuring m^, at low energy only 
determines 9 of the 18 seesaw parameters. We will see below that leptogenesis and lepton 
flavor violation are related to different combinations of the 18 seesaw parameters and can 
provide different information to determine the seesaw parameters from the z/-oscillation and 
LEV observations. 

We can rewrite the seesaw formula Eq. ()19j) in another form 

U.J^Ju^J^Y = -Vm%X^^=\Vm%X^^=] , (21) 

from which rriN can be solved in terms of the left- and right-handed neutrino masses. 



m 



N 



Vymi0^jMf,X* , (22) 



where O is an arbitrary orthogonal 3x3 matrix |3 71] and V = UlV . In the above equation 
we have absorbed all the 6 Majorana CP phases in the diagonal eigenvalue matrices: two 
low energy Majorana phases, r/, ^, are absorbed by J ml and the four high energy Majorana 



phases, a, /3, p, uj, are absorbed by and y Mf^. We will use this equation repeatedly in 
the following discussions. 



B. The reheating temperature 

The lightest sneutrino A'^i begins to oscillate when the Hubble expansion rate H ~ Mi 
and decays at H ~ P^v ■ The Universe is then reheated by the relativistic decay products. 
The reheating temperature is approximately determined by 

T^H-(^^j ^r^Mp , (23) 

where g^^ is the number of the effective relativistic degrees of freedom in the reheated Uni- 
verse, Mp = I/v^SttG^ ~ 2.4 X lO^^GeV^ is the Planck scale, and 

= ^(4r^)iiMi , (24) 
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is the width of the hghtest sneutrino A^i, if it couples to other matter only through the 
Yukawa coupling in Eq. Taking Mi ^ 1.7 x lO^^GeV and Trh ~ lO^^GeV, we get 

(Y^Y/v)ii should be as small as 0(10"^°). 

The reheating temperature (as well as leptogenesis) is related to the RH mixing of Yn 
and put strong constraints on this mixing matrix. Using Eq. (fT7|) . we have 



(25) 



The elements \Xii\ can be parametrized by two mixing angles, 61^2- We then get 



(4y^)n = cjclY,' + clslYi + sjYi ^ Y^ + + s\Y,' ^ 10" 
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(26) 



with Ci = cosdi, Si = sin6'j. In the later discussion we will see that Y2 is (9(0.1) and Y3 is 
(9(1). Then we have 

< 10-^° , si < 10-« , sj < 10-^° . 



(27) 



Since 6'i 2 are extremely small, X can be given in a quite simple form as 



X 



1 S2 Si 

-ic3S2 + S3Sl) C3 S3 
V S3S2 - C3S1 -S3 C3 J 



(28) 



where Si = Sie*^. 

Using Eqs. ((22) and (jH)), we have 



(mjvmAr)ii 



{v sin p)^ ^ 4 X 10-^ 



T, 



RH 



Mim^AOnl'^ + UOJGeV^lOuf + 850GeV^\di3\'^ 



(29) 



where we have assumed sin /5 ^ 1 for large tan /5, and m,. 



^"^L ~ 7.1 X lO^^el^ and 



Am; 



atm 



0.05e\^. From the above equation we can see that Ou and O13 have to 



be negligibly small. We will set these two elements zero and write O as 

(±1 \ 

6 = 



V 



c s 
— s c 



(30) 



where c = cosOt, s = sin 6't with 6t being an arbitrary complex angle. (It should be 
noted that O12 and O13 can not be exactly zero, since if they are zero the first-generation 
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right-handed (s)neutrino decouples from the other two generations and no lepton number 
asymmetry can be induced when it decays. However, the tiny mixing has no effect on lepton 
flavor violation and we can ignore them safely when discussing LFV.) 
From Eq. (j29|) we can estimate that 

f 2. X 10-i°eV, for Trh = lO^^GeV , ^ ^ 

rrii,^ ^ < (31) 
[ 2. X for Trh = lO^GeV . 

This estimation is correct when the last two terms are much smaller than the first one in 
the second line of Eq. (j^^ . or, equivalently, the Yi term dominants the others in Eq. ()26p. 
In the following discussion for leptogenesis we will see that this is a quite natural situation. 



C. Leptogenesis 



Since the reheating temperature, Trh, is far below the lightest RH (s)neutrino mass. Mi, 
leptogenesis arises dominantly from direct cold sneutrino decays, with negligible thermal 
wash-out effects. In this case, the baryon asymmetry is given byjssl 

_ ub 3 Trh . . 

where a = —8/23 is the ratio of baryon to lepton asymmetry balanced by the "sphaleron" 
process. In order to produce the observed baryon asymmetry in the Universe, Yb ~ 10^^°, 
we require the sneutrino decay asymmetry ei ~ —10"^. The asymmetry ei is given by 

^1 - -4z7TA^ E Ini liYlY^uT^ . (33) 



levr (r^r^)i 



^ 1 1=2,3 



Mi 

Using the expression for X in Eq. ()28|) and the large hierarchy between Yi and Y2, Y3 we 

get 

{Y^Y^)^2 ^ {s2Yi + slssY.y^ , 

{yIYn)i3 ^ {-S2S3Yi + slY^)e''' . (34) 

We will discuss two simple cases to illustrate some quantitative features of the seesaw 
parameters required by leptogenesis. We will see that, in Eq. fOF)|) . the Y2 and Y3 terms 
should be smaller than the Yi term in order to produce the lepton number asymmetry at 
the correct order. 
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Case I, siY^ < 52^2^ 

In this case the expression for ei is simphfied as 



167r (4Fjv)ii 

- -10-^ -TT^^- sin 2a. (35) 

When deriving the second hne we have assumed that a and (3 are of the same order 
and ^ < ^ ~ 1^2 ~ 0.1 and sjYi < ^2^22. if the CP phases are of order 1, s^Y^/Y^ 
should be at the order of about 10^^. Actually, this case corresponds to the maximal 
asymmetry given by ^ -|_:^^£iv^^£°i ^ 10~^ j^^. In this case, the CP phase or, 

slYi/{YlYN)ii, has to be at the order of C(10-2). 

Case II, si ~ S2 

In this case we can simplify the expression for ei as 



3 1 



{s,s,Yirsm2a'^^+is,Yirsm2(3'^^ 



167r (4r^)n 

10"^ ■ v2 ' \v2 ■ (sin 2a' + sin 2/3') , (36) 

where we have used the fact that S3 ~ -1/M2/M3 if 9t is of order 1, sfY^ ^ SIY2 and 



a' = a — 6, P' = P — 6. Similar to Case I, we get that s^Y^/Y^ should be at the 
order of 10~^ if the CP phases are of order 1. In this case, the maximal asymmetry is 

\^max \ ^ 3 A/iy/Amg^^ 3 
I J- I Idtt 

Certainly, it is possible that the contributions to (yIjYn)ii from I2 and F3 in Eq. (pUj) 
are of the same order. In this case we also expect that these values be correct as an estimate 
of the order of magnitude, i.e., slY^ ~ s\Y2 ^ Y^. This analysis justifies our guess in the 
last subsection that the Yi term gives the dominant contribution in the process of reheating 
the Universe. Conversely, if the Y2 or Y3 term gives dominant contribution, the CP phases 
have to be fine tuned to the order of 10"^ and 10"'^ respectively, in order not to create too 
much lepton number asymmetry and m^j^ in Eq. (j^T|) will be even smaller. 



D. Lepton flavor violation and muon anomalous magnetic moment 



We have shown that leptogenesis is associated with the high energy mixing angles and 
CP phases in the unitary matrix X. Generally, leptogenesis has no direct relation with the 
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low energy neutrino phenomena. However, another interesting phenomena — the charged 
lepton flavor violating decays — predicted by this sneutrino inflaton model, can provide 
constraints on the seesaw model's parameter space. The muon anomalous magnetic moment 
is also considered to constrain the SUSY parameters. 

In a supersymmetric model, the present experimental limits on the LFV processes has put 
very strong constraints on the soft supersymmetry breaking parameters, with the strongest 
constraints coming from the process /i — > 67 (BR(/i 67) < 1.2 x 10~^^j4^). It is a 
usual practice to assume universal soft SUSY breaking parameters mg, Tni/2 and Aq at 
the SUSY breaking scale ( We take it the GUT scale here) to suppress the LFV effects. 
However, since there are LFV interactions in the seesaw models, the lepton flavor violating 
off-diagonal elements of (m|)jj, the slepton doublet soft mass matrix, and {Ae)ij, the lepton 
soft trilinear couplings, can be induced when running the renormalization group equations 
(RGBs) for m| and between Mqut and Mr. 

The off-diagonal elements of {Tn?i)ij and {Ae)ij can be approximately given by 

(K),,. ^ ^(^^4)..(3 + a^Wo log ^ , (37) 

(Me),,. ^ ^(r^4),,amolog^^ , (38) 
where Aq = amo is the universal trilinear coupling at Mqut- Using Eq. (fTTj) we have 

{YnYI\, = {V{Y^fV% ^ Vi2V;,Yi + V,sV;,Y,' , for Mgut > Q > Ms (39) 

k=l,2 

- E VuV*^YiYUSim - X3iX;j , for Ms > Q > M2 . (40) 

l,m=2,3 

The numerical result shows that, since the mixing angles in X are all small, the LFV effects 
are only sensitive to the left-handed mixing matrix V, while leptogenesis only relies on 
the right-handed mixing matrix X. Thus, there are no direct relation between the two 
phenomena in principle. 

We have solved the full coupled RGEs numerically from the GUT scale to Mz scale. 
At the energy scales below M2 we solve the RGEs for MSSM and below Msusy the RGEs 
return to those of the SM. 

In principle, only 9 of the 18 seesaw parameters are determined in our model, i.e., m,^., 
Mj and 3 low energy neutrino mixing angles. In order to predict the branching ratio of the 
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LFV decays, we have to explore a 9-dimensional parameter space of the unknown variables. 
However, from our previous discussions, we know that the relevant seesaw parameters to 
LFV are reduced to only 4 in this model, which can be chosen as 1 complex angle Ot-, and 
2 CP phases. We can explicitly write the 'reduced' seesaw formula for the 2nd and 3rd 
generations as 

/ 



\ 




(41) 



where both V and X are 2x2 real orthogonal matrices, determined by diagonalizing the 
matrix M.. Here, we adopt the running values of niy^ and at the scale of IQ^'^GeV 41 1. 
Once 1^2,3 and V = U^V are determined, can we calculate the LFV branching ratios, BR(/j — »• 

hi)- 

The relevant parameters to investigate BR(/j Ij'y) and 6a^ include the mSUGRA 
parameters: mo, mi/2, Aq, tan/5, sgn(/i) and the seesaw parameters: 9t and 0j. Since 
BR(Zj — > Ijj) and nearly scale with tan^ (3 and tan jS respectively, we take tan /5 = 10 as 
a representative value. We fix = through our calculation since it has small influence 
on the numerical results. The Higgsino mass parameter /x > is assumed, motivated 
by the ~ 2 anomaly. As for the seesaw parameters, we take Am^^^ = 5 x 10~^eV^^, 



Am' 



atm 



2.5 X 10"3e\/2, and tan^ 6112 = 0.42, sin^ 26'23 = 1, < 6113 < 0.2 from the neutrino 



oscillation experiments. We fix Mi = 1.7 x lO^GeV, M2 = lO^^GeV" and 4 x W^GeV < 
M3 < 1 X lO^^GeV for RH heavy Majorana neutrinos. 

In Fig. 13 we plot BR(/j Ij'j) and Sa^ as functions of mi/2 and mo for 6t = 7r/4 and 
(pi = 0. From this figure we can see that the process fi —>■ ej gives very strong constraint on 
the SUSY parameter space: only with large mi/2 and relatively small mo can its branching 
ratio be below the present experimental limit, 1.2 x 10^^^. For the following discussions, we 
will fix (mi/2,'^o) = (800, 250)Gel^. Since the muon anomalous magnetic moment, 6a^, is 



42l |. it is also fixed at about 2.7 x 10 which 



nearly independent of the seesaw parameters 
will be omitted in the other figures. 

Taking determinant on Eq. 1)411) we know that the product of 1^2,3 is fixed by the left- 
and right-handed Majorana neutrino masses. The ratio of the two Yukawa couplings is 
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FIG. 7: BR(/j Ijj) and as a function of mi/2 and mo in the left and right panels respectively, 
tan/? = 10, Aq = and /i > are fixed. We take mo = 250GeV for the left panel and = 
SOOGeV for the right panel. The seesaw parameters are taken as 9t = 7r/4, (pi = 0, ^13 = 0.05 and 
Ms = 1 X lO^^GeV. 
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Re(6»T) Re(6»T) 

FIG. 9: BR(/i ^ 67) (Left) and BR(r ^ /i7) (Right) as function of ReOr for Im6'r = 0, 0.5, 1, 1.5. 
We fix 6*13 = 0.05, (pi = and M3 = 1 x lO^^GeV. 

determined by djn. In Fig. |Hlwe show Y3 as function of KedT and ImOx- Both the real and 
imaginary part influence the ratio between Y2 and 1^3. Since 1^3 increases almost linearly 
with Im6'7-, we expect BR(/j Ijj) also increase with Im6'j'. 

In Fig. ini we plot BR(/i — > 67) and BR(r — > /i7) as function of Re6'r on the left and right 
panels respectively. For Im^j- = 0.5, BR(yU 67) has been greater than the experimental 
limit. 

In Fig. BR(/j — ^ Ij'y) is drawn as function of ^13. We can see BR(/i — 67) is very 
sensitive to ^13, while BR(r —>■ fi'-f) is insensitive to ^13. The behavior in this figure is 
understood if we notice that the flavor mixing between the first and the second generations 
is nearly proportional to ^^13^*31^3^, where V13 = (?7i/)i2V^3 + (f^i/)i3^33- The two terms are 
added constructively or destructively, depending on the sign of 6'^. When we set 612 = 0, 
the branching ratio of /i — 67 increases rapidly with ^13, independent of the value of 6t- 

In Fig. [m we plot BR(/j Ij^y) as function of 0i, which determines the relative phase 
between f/j, and V. The behavior in the figure is easy to understand. We also examined 
that BR(/i Ij^) is indeed independent of 02, as we expected. 
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^13 

FIG. 10: BR(/i ^ Ijj) as function of 613. We &x cp, = and M3 = 1 x lO^^GeF. 

Finally, we plot BR(/j —>■ Ij'-f) as function of M3. BR(/j —>■ Ijj) increases with M3 at 
first, because it makes 12,3 larger. However, when M3 is as large as lO^^GeV, which is too 
close to Mquti the integration distance log ^^^^^^ becomes too small and the branching ratio 
decreases. Although below M3 LFV is still produced, see Eq. fjlUjl . the effects are small, 
since the contribution from Y2 is small, due to Y2 « Y.^. The 13 coupling contributes 
to the LFV below M3 through the mixing, Y3^|X23p, which is also small due to the small 
mixing element. 

We have omitted BR(r 67) in all the figures because the predicted branching ratio is 
much smaller than the present experimental limit. 

E. Neutrinoless double beta decay 

From the above discussion we know that it is impossible to produce degenerate solution 
for the left-handed neutrino masses in this model. It is easy to estimate that < m >ee= 
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FIG. 11: BR(/i ^ Ij-i) as function of (pi. We fix 6*13 = 0.05, 02 = and M3 = 1 x lO^^GeF 



(2 ~ 4) X 10 ^eV, depending on the value of f/es- So, in this sneutrino-inflaton model, it is 
hard to account for the neutrinoless double beta decay experimental signalj28|. 



IV. SUMMARY AND DISCUSSIONS 



We have considered a double-sneutrino inflation model within the minimal supersymmet- 
ric seesaw model. With the mass ratio 6 ^ r ^ 10 and the lighter sneutrino Mi ~ 1.7 x 10^^ 
GeV, the model predicts a suppressed primordial scalar spectrum around the largest scales 
which is favored at > 1.5a. The predicted CMB TT quadrupole is much better suppressed 
than the single sneutrino model and the preference level by the WMAP first year data is 
about 1.5a. Double quartic inflation can also work very well in light of WMAP observations. 

We then have studied the phenomenological implications of this model. The seesaw 
parameters are constrained by both particle physics and cosmological observations. The 
strongest constraint comes from the required reheating temperature by the gravitino prob- 
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FIG. 12: BR{li Ijj) as function of M3 for 0t = 7r/4 + 0.5i and Or = it/4:. We fix 6*13 = 0.05 
and (pi = 0. 

lem. To some extend, fine tunning is needed to satisfy this constraint, which means that 
the right-handed mixing angles 6i and 62 are much smaller than the mass hierarchy of the 
right-handed neutrinos. Further, the mass of the lightest left-handed neutrino should be at 
the order of 10~^^eV, much smaller than the other two light neutrinos. 

Leptogenesis arises from the decays of the cold infiaton — the lightest sneutrino. It is 
easy to account for the observed quantity of the baryon number asymmetry in the Universe 
by adjusting the seesaw parameters. 

This model gives definite predictions on the lepton flavor violating decay rates. In most 
parameter space, the branching ratio of /x — > 67 is near or exceeds the present experimental 
limit. However, the branching ratio of r — > yU7 is at the order of about 10~^° — 10~^, which 
is far below the current experimental limit. Furthermore, in the appropriate range of SUSY 
parameter space where LFV constraints are satisfied, the SUSY can only enhance the muon 
anomalous magnetic moment at the amount of (2 ~ 3) x 10~^°. 
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This model can not predict a degenerate light neutrino spectrum. The observed signal 
of neutrinoless double beta decay, if finally verified, can not be explained by the effective 
Majorana neutrino mass in this model. 

Acknowledgments 

We thank Z. Z. Xing for helpful discussions. We acknowledge the using of CMBFAST 
program This work is supported by the National Natural Science Foundation of 

China under the grant No. 10105004, 19925523, 10047004 and also by the Ministry of 
Science and Technology of China under grant No. NKBRSF G19990754. 



[1] A. Guth, Phys. Rev. D 23, 347 (1981); A. Linde, Phys. Lett. B 108, 389 (1982). A. Albrecht 

and P. J. Steinhardt, Phys. Rev. Lett. 48, 1220 (1982). 
[2] H. Murayama, H. Suzuki, T. Yanagida, and J. Yokoyama, Phys. Rev. Lett. 70, 1912 (1993); 

H. Murayama, H. Suzuki, T. Yanagida, and J. Yokoyama, Phys. Rev. D 50, 2356 (1994). 
[3] J. Ellis, M. Raidal and T. Yanagida, arXiv: |hep-ph/0303242 



[4] Y. Fukuda et al. [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81 (1998) 
1562; Q. R. Ahmad et al. [SNO Collaboration], Phys. Rev. Lett. 89 (2002) 011301, 



arXiv:nucl-ex/0204008 K. Eguchi et al. [KamLAND Collaboration], Phys. Rev. Lett. 90 



(2003) 021802, arXiv:hep-ex/0212021 



[5] M. Gell-Mann, P. Ramond and R. Slansky, Proceedings of the Supergravity Stony Brook 
Workshop, New York, 1979, eds. P. Van Nieuwenhuizen and D. Freedman (North-Holland, 
Amsterdam); T. Yanagida, Proceedings of the Workshop on Unified Theories and Baryon 
Number in the Universe, Tsukuba, Japan 1979 (eds. A. Sawada and A. Sugamoto, KEK 
Report No. 79-18, Tsukuba); R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 
(1980); S. L. Glashow, Caraese lectures, (1979). 

[6] M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45; P. Langacker, R. D. Peccei 
and T. Yanagida, Mod. Phys. Lett. A 1, (1986) 541; R. N. Mohapatra and X. Zhang, Phys. 
Rev. D 46, (1992) 5331; H. Murayama and T. Yanagida, Phys. Lett. B 322 (1994) 349 
arXiv:hep-ph/9310297| K. Hamaguchi, H. Murayama and T. Yanagida, Phys. Rev. D 65 



25 



(2002) 043512 |arxrv:hep-ph/01Q9030| T. Moroi and H. Murayama, Phys. Lett. B 553 (2003) 



126 arXiv:hep-ph/0211019 



[7] C. L. Bennett et a^., |astro-ph/0302207, 
[8] G. Hinshaw et aL, |astro-ph/0302217[ 
[9] C. L. Bennett et ai, Astrophys. J. 464, LI (1996). 
[10] H. V. Peiris et aZ., |astro-ph/0302225( 

[11] W. J. Percival et al, Mon. Not. Roy. Astr. Soc. 327, 1297 (2001). 
[12] P. Mukherjee and Y. Wang, astro-ph/03032l'Tj 

[13] S. L. Bridle, A. M. Lewis, J. Weller, and G. Efstathiou, |astro-ph/Q30'2306| 
[14] Y. Wang, D. N.Spergel, and M. A.Strauss, Astrophys.J. 510,20 1999; Y. Wang and G. Math- 
ews, Astrophys. J. 573,1 2002; P. Mukherjee and Y. Wang, astro-ph/0301058 ; P. Mukherjee 
and Y. Wang, astro-ph/0301562l| 



[15] D. N. Spergel et al, |astro-ph/0302209 

[16] J. M. Cline, P. Grotty and J. Lesgourgues, |astro-ph/0304558| G. Efstathiou, [astro-ph/030643H 
A. d. Oliveira-Costa, M. Tegmark, M. Zaldarriaga and A. Hamilton, astro-ph/0307282j A. 



Niarchou, A. H. Jaffe and L. Pogosian, astro-ph/030846lj 



[17] Y.-P. Jing and L.-Z. Fang, Phys. Rev. Lett. 73,1882 (1994); J. Yokoyama, Phys. Rev. 
D59, 107303 (1999); S. DeDeo, R. R. Caldwell and P. J. Steinhardt, Phys.Rev.D 67, 
103509 (2003) ; J. Uzan, A. Riazuelo, R. Lehoucq and J. Weeks, |astro-ph/0303580| ; 
G. Efstathiou, Mon.Not.Roy.Astron.Soc. 343, L95 (2003); G. R. Contaldi, M. Peloso, 
L. Kofman, and A. Linde, JCAP 0307, 002 (2003); E. Gaztanaga, J. Wagg, T. Mul- 



hashi, |hep-ph/0305319} M. Bastero-Gil et ai, in Ref. 



astro-ph/030'^ 


117 


'8 


ai, in Ref. 


21 


1; 



2l|]; A. Lasenby and C. Doran, 



astro-ph/03073 1 1 ( S. Tsujikawa, R. Maartens and R. Brandenberger, astro- ph/0308169| T. 



Moroi and T. Takahashi |astro-ph /0308208 ; Q. Huang and M. Li, |astro-p h/0308458, 
[18] B. Feng and X. Zhang, Phys.Lett.B 570, 145 (2003) 
[19] B. Feng, M. Li, R.-J. Zhang, and X. Zhang, |astro-ph/0302479 
[20] X. Wang et a?., astro-ph/0209242] 



[21] J. E. Lidsey and R. Tavakol, astro-ph /03041 13| M. Kawasaki, M. Yamaguchi and J. Yokoyama, 
Phys.Rev. D 68, 023508 (2003); Q. G. Huang and M. Li, JHEP 0306, 014 (2003); D. J. 
Chung, G. Shiu and M. Trodden, |astro-ph/0305l"93l K.-I. Izawa, |hep^^0305286 ; M. Bastero- 



26 



Gil, K. Freese and L. Mersini-Houghton, hep-ph / 03Q6289| M. Yamaguchi and J. Yokoyama, 
pIep^ph/03073731 

[22] V. Barger, H. Lee, and D. Marfatia, Phys.Lett.B 565, 33 (2003); W. H.Kinney, E. W.Kolb, A. 

Melchiorri and A. Riotto, piep'-ph/0305130l S. M.Leach and A. R.Liddle, jastro-ph/0306"305| 
[23] D. Polarski and A. A. Starobinsky, Nucl.Phys.B 385, 623 (1992); D. Polarski, Phys.Rev.D 49, 

6319 (1994); D. Polarski and A. A. Starobinsky, Phys.Lett.B 356, 196 (1995); J. Lesgourgues 

and D. Polarski, Phys.Rev.D 56, 6425 (1997). 
[24] M. H. Kesden, M. Kamionkowski and A. Cooray, ,astro-ph /0306597^ 
[25] M. Kamionkowski, in private communications. 

[26] J. M. Diego, P. Mazzotta and J. Silk, |astro-ph/0309181[ O. Dore, G. P. Holder and A. Loeb, 



[astro-ph /0309281; P. G. Castro, M. Douspis and P. G. Ferreira, |astro-ph /0309320| 
[27] J. R. Ellis, J. E. Kim and D. V. Nanopoulos, Phys. Lett. B 145 (1984) 181; J. R. Ellis, 
D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 259 (1985) 175; J. R. Ellis, D. V. Nanopou- 
los, K. A. Olive and S. J. Rey, Astropart. Phys. 4 (1996) 371; M. Kawasaki and T. Moroi, 
Prog. Theor. Phys. 93 (1995) 879; T. Moroi, Ph.D. thesis, |arXiv:hep-ph/9503210| M. Bolz, 

A. Brandenburg and W. Buchmiiller, Nucl. Phys. B 606 (2001) 518; R. Cyburt, J. R. Ellis, 

B. D. Fields and K. A. Olive, |astro-ph/0211258 



[28 
[29 
[30 

[31 
[32; 
[33 
[34; 
[35; 
[36 
[37 
[38 



H. V. Klapdor-Kleingrothaus et al, Mod. Phys. Lett. A 16, 2409 (2002). 

C. Gordon, D. Wands, B.A. Bassett, and R. Maartens, Phys. Rev. D63, 023506 (2001). 

S. M. Leach and A. R. Liddle, Phys. Rev. D 63, 043508 (2001); S. M. Leach, M. Sasaki, 

D. Wands and A. R. Liddle,Phys. Rev. D 64, 023512 (2001). 

D. H. Huang, W. B. Lin and X. M. Zhang Phys. Rev. D 62, 087302 (2000). 



B. Feng, X. Gong and X. Wang, astro-ph/0301111 



D. H. Lyth and A. Riotto, Phys.Rept. 314, 1 (1999). 



L. Verde et al, arXiv: |astro- ph/0302218| 

A. R. Liddle and S. M. Leach, arXiv: astro-ph/0305263l 

Z. Maki, M. Nakagawa, and S. Sakata, Prog. Theor. Phys. 28, 870 (1962). 

J. A. Casas and A. Ibarra, Nucl. Phys. B 618, 171 (2001). 

K. Hamaguchi, H. Murayama, T. Yanagida, Phys. Rev. D 65, 043512 (2002), arXiv: 



|hep-p h/0109030| 

[39] W. Buchmiiller, P. Di Bari and M. Pliimacher, Nucl. Phys. B 643, 367 (2002). 



27 



[40] M. L. Brooks et ai, MEGA Collaboration, Phys. Rev. Lett. 83 1521 (1999); S. Ahmed et al, 
CLEO Collaboration, Phys. Rev. D 61 071101 (2000); U. Bellgardt et al, Nucl. Phys. B229 
1 (1988). 

[41] S. Antusch, J. Kersten, M. Lindner, M. Ratz, arXiv: |hep-ph/0305273| 
[42] X. J. Bi, Y. P. Kuang and Y. H. An, arXiv: hep-ph/021 11421 
[43] U. Seljak and M. Zaldarriaga, Astrophys. J. 469, 437 (1996). 
[44] |http://cmbfast.org/ . 



28 



